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Herein,  spherical  nanoindentation  (NI)  was  used  to  investigate  the  room  temperature  deformation 
behavior  of  C-plane  LiTa03  single  crystals  loaded  along  the  [0001]  direction  as  a  function  of  ion 
irradiation.  When  the  NI  load-displacement  curves  of  3  different  nanoindenter  radii  (1.4  /(in,  5  /mi, 
and  21  /mi)  were  converted  to  NI  stress-strain  curves,  good  agreement  between  them  was  found. 

The  surface  first  deforms  elastically  -  with  a  Young’s  modulus  of  205  ±  5  GPa,  calculated  from 
the  stiffness  versus  contact  radii  curves  and  207  ±  3  GPa  measured  using  a  Berkovich  tip  -  and 
then  plastically  deforms  at  «6  GPa.  Repeated  loading  into  the  same  location  results  in  large, 
reproducible,  fully  reversible,  nested  hysteresis  loops  attributed  to  the  formation  of  incipient  kink 
bands  (IKBs).  The  latter  are  coaxial  fully  reversible  dislocation  loops  that  spontaneously  shrink 
when  the  load  is  removed.  The  IKBs  most  probably  nucleate  within  the  ( 1  Ol 2)  twins  that  form  near 
the  surface.  The  sharper  radii  resulted  in  twin  nucleation  at  lower  stresses.  The  changes  in  the 
reversible  loops’  shape  and  areas  can  be  related  to  the  width  of  the  twins  that  form.  The  latter  were 
proportional  to  the  nanoindenter  tip  radii  and  confirmed  by  scanning  electron  microscopy  and  by 
the  fact  that  larger  threshold  stresses  were  needed  for  IKB  nucleation  with  the  smaller  tip  sizes.  No 
effect  of  irradiation  was  observed  on  the  NI  response,  presumably  because  of  the  mildness  of  the 
irradiation  damage.  ©2011  American  Institute  of  Physics,  [doi:  10. 1063/1. 3608 158] 


I.  INTRODUCTION 

Lithium  tantalate,  LiTaOs,  is  one  of  the  most  attractive 
materials  for  nonlinear  integrated  optics.  The  properties  of  fer¬ 
roelectric  LiTa03  are  similar  to  those  of  lithium  niobate, 
LiNb03,  reviewed  in  detail  elsewhere,1’2  with  excellent 
nonlinear  optical  properties,  as  well  as  higher  thresholds  to 
photorefractive  damage  than  LiNb03  and  a  lower  Curie  tem¬ 
perature.’  The  importance  of  LiTa03  is  emerging  in  applica¬ 
tions,  such  as  second-harmonic  generation  and  optical 
parametric  amplification/oscillation,  because  it  can  be  proc¬ 
essed  to  make  waveguides  and  to  engineer  the  nonlinearity 
through  quasi-phase  matching  via  periodic  poling,4  offering 
an  important  alternative  to  LiNb03  in  devices  for  signal  proc¬ 
essing  via  quadratic  cascading.5  Despite  these  potential  appli¬ 
cations,  little  work  has  been  carried  out  on  its  mechanical 
behavior  especially  at  room  temperature.  Doukhan  et  al.6 
described  different  lattice  defects  in  LiTa03  and  concluded, 
not  surprisingly,  that  it  had  the  same  twinning  system  as 
LiNb03. 

Recently,  we  showed  that  the  vast  majority  of  plastically 
anisotropic  solids  with  c/a  ratios  >1.4  can  be  classified  as 
kinking  nonlinear  elastic  (KNE)  solids. 7-10  A  sufficient  con¬ 
dition  for  a  solid  to  be  KNE  is  plastic  anisotropy.  The  signa¬ 
ture  of  these  solids  is  the  formation  of  fully  reversible, 
reproducible  stress-strain  loops  during  cyclic  loading.  The 
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full  reversibility  of  these  loops  is  believed  to  be  caused  by 
incipient  kink  bands  (IKBs)  that  are  comprised  of  multiple, 
coaxial,  parallel  dislocation  loops  (Fig.  1(a)),  which  remain 
extended  only  if  the  load  is  applied;  when  the  load  is 
removed,  they  shrink  or  are  annihilated  altogether. 1 1  At  high 
stresses,  such  as  under  a  nanoindenter,  the  IKBs  sunder  and 
devolve,  first  into  mobile  dislocation  walls  (MDW’s)  and 
ultimately  into  kink  boundaries  (KBs),  which  are  irreversible 
(Fig.  l(b)).7’10'12~is  On  reloading  to  the  same  stress,  IKBs 
nucleate  in  the  newly  created  microdomains.  The  to-and-fro 
motion  of  the  IKB  dislocations,  in  turn,  results  in  hysteresis 
and  the  dissipation  of  energy.  Using  this  approach,  we 
showed  that  graphite,15  mica,7'17  ZnO  (C-orientation),12  sap¬ 
phire,13  and  GaN,14  among  many  others,  are  KNE  solids. 

Following  the  pioneering  work  of  Herbert  et  al.19  and 
Oliver  and  Pharr,  Suganuma,  and  Bushby,  “  we  devel¬ 
oped  a  technique  for  converting  spherical  NI  load/displace¬ 
ment  curves  to  NI  stress/strain  curves.23’24  According  to  our 
technique,  NI  stress  and  strain  are  defined  as  P/7ia2  and  a/R, 
respectively,  where  P,  a,  and  R  are  load,  contact  radius,  and 
tip  radius.  By  plotting  the  NI  stresses  and  strains,  more  infor¬ 
mation  can  be  gleaned  from  the  indentation  results.  We 
applied  this  method  to  better  understand  the  deformation 
behavior  of  a  number  of  oxide  and  nitride  single  crystals, 
such  as  sapphire,13  ZnO,12  LiNb03,25  mica,7'17  GaN,14 
BaTKU,26  and,  more  recently,  polycrystalline  Y2O3.27 

Most  relevant  to  this  paper  is  our  recent  work  on  the 
nanoindentation  of  C-plane  LiNb03  single  crystals  loaded 
along  [0001]. "5  In  that  paper,  we  reported  an  elastic  modulus 
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(B)  Mobile  Dislocation  Walls 


FIG.  1 .  (Color  online)  Schematic  of  (a)  dislocation  loops  comprising  an  IKB; 
(b)  schematic  of  what  could  be  occurring  below  the  indented  surface.  The 
emission  of  mobile  dislocation  walls  that,  in  turn,  form  kink  boundaries,  are 
shown.  Upon  re-loading,  the  IKBs  form  within  the  kink  boundaries  or  twins 
formed  during  the  pop-in  events,  (c)  Typical  stress-strain  curve  for  a  KNE 
solid  obtained  under  spherical  nanoindentations,  showing  the  definition  of  non¬ 
linear  strain,  e^L,  and  the  energy  dissipated  per  unit  volume  per  cycle,  Wd. 


of  186  GPa  and  a  Vickers  microhardness  of  about  4.5  GPa. 
Cyclic  loading  resulted  in  the  signature  of  KNE  solids,  viz., 
large,  fully  reversible,  reproducible,  hysteretic  stress-strain 
loops.  As  far  as  we  are  aware,  these  remain  the  largest  ever 
reported  for  crystalline  solids. 

Before  discussing  the  work  carried  out  herein,  it  is  im¬ 
portant  to  summarize  our  IKB-based  model  that  is,  in  turn, 
based  on  early  work  by  Frank  and  Stroll  (F  &  S).11  The  fol¬ 
lowing  is  a  simplified  version.  F  &  S  considered  an  elliptic 
kink  band  (KB)  with  length,  2a,  and  width,  2/1,  such  that 
a  /I  (Fig.  1(a))  and  showed  that  the  remote  shear  stress,  t, 
needed  to  render  such  a  subcritical  KB  unstable  is  given  by 


M  V  7t/2a  wyc 


(1) 


where  tc  and  <r,  are  the  remote  critical  shear  and  axial 
stresses,  resectively.  M  is  the  Taylor  factor  relating  the  shear 
stress  at  the  grain  level  to  the  applied  stress.  The  maximum 
value  of  M  is  2.  G,  b,  and  w  are,  respectively,  the  shear  mod¬ 
ulus,  Burgers  vector,  and  a  term  related  to  the  dislocation 
core  width. 1 1  If  one  assumes  that  the  local  stress  needed  to 
nucleate  an  1KB  is  «  G/35,  then  at  0.05  rad.,  the  critical 
shear  angle,  yc,  is  small.11’28 

As  a  first  approximation,  each  dislocation  loop 
(Fig.  1(a))  can  be  assumed  to  be  comprised  of  two  edge  and 
two  screw  dislocation  segments  with  lengths  2/jx  and  2/>y, 
respectively.  It  is  also  assumed  that,  when  cr  >  cr,,  the  IKBs 
grow  by  increasing  their  width,  2/1,  according  to 


2/1, 


2«(!  ~  v)  g  9/? 
Gyc  M'  ^ 


2a  cr 
GycM 


(2) 


for  the  edge  and  screw  components,  respectively.  It  follows 
that,  for  cr  >  crt,  the  IKBs  grow  and  the  IKB-induced  axial 
strain  resulting  from  their  growth  is  assumed  to  be  given  by29 


£ikb  — 


AVNkyc  __  Nkyc4nct{PJy  - 
k\ 

4n(l  —  v)Nk  a3 


3kiG2ycM 2 


3ki 

(o-2  -  ff/)  =  mj(f72  -  cr2),  (3) 


where  m,  is  the  coefficient  before  the  term  in  brackets  in  the 
fourth  term,  Nk  is  the  number  of  IKBs  per  unit  volume,  and 
AV  is  the  change  in  the  volume  kinked  as  the  IKBs  grow 
from  a  size  at  cr,  to  their  size  at  cr.  The  factor  k,  relates  the 
volumetric  strain  due  to  the  IKBs  to  the  axial  strain  along  the 
loading  direction.  Experimentally,  k,  varies  from  1  to  2.  For 
example,  in  polycrystalline  Mg,  it  is  closer  to  l,30  while  in 
Co,  it  is  closer  to  2. 31  Reed-Hill  et  at ?2  also  assumed  k,  =2 
when  they  modeled  twins  in  Zr.  Herein,  we  assumed  k,  =  2; 
the  implications  and  ramifications  of  this  assumption  are  dis¬ 
cussed  below.  Once  m,  is  determined  experimentally,  if  2a 
can  be  estimated,  Nk  can  be  calculated.  Note  that  Nka3  is  of 
the  order  of  unity. 

If  Q  is  the  energy  dissipated  by  a  dislocation  line  sweep¬ 
ing  a  unit  area,  then  the  energy  dissipated  per  cycle  per  unit 
volume,  Wd,  can  be  expressed  as8'29 


Wd 


m2(cr2  -  cr2), 


PM 


3k\—m\(o2  -  cr2) 

(4) 


in  which  D  is  the  distance  between  dislocation  loops  along 
2a  (Fig.  la).  It  follows  that  Q/b  should  be  proportional,  if  not 
equal,  to  the  critical  resolved  shear  stress  (CRSS)  of  an  1KB 
dislocation  loop.  Combining  Eqs.  (3)  and  (4)  yields 

Q  m2 

Wd  =  3k\—£iKB= — e ikb-  (5) 

b  m\ 

Figure  1(b)  shows  a  schematic  of  how  £nl  and  Wd  are  esti¬ 
mated  from  the  NI  stress-strain  curves. 

Assuming  the  IKBs  to  be  cylinders  with  radii  /)av,  then 
the  reversible  dislocation  density,  prev,  due  to  the  IKBs  is 
given  by 
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<•'«  =  — — = — b — ■  (6> 

where  /Jav  is  the  average  of  j8x  and  f>y. 

Since  IKBs  are  comprised  of  dislocation  loops  and  NI  is 
sensitive  to  near  surface  properties,  it  is  reasonable  to 
assume  that  any  process  that  induces  near  surface  defects 
could  affect  the  NI  response  of  KNE  solids.  This  is  particu¬ 
larly  true  since  we  recently  showed  that  the  average  pop-in 
stresses  in  defective  micas  were  significantly  lower  than 
those  in  less  defective  ones.17  At  this  time,  it  is  fairly  well 
established  that  irradiation  of  solids  with  light  ions,  such  as 
helium  (He),  can  create  a  large  number  of  defects,  including 
vacancies  and  interstitials,  voids,  dislocation  loops,  and  He 
bubbles,  etc.3’  34  On  the  other  hand,  there  are  some  parame¬ 
ters,  such  as  the  material’s  compositional  complexity,  that 
can  act  to  suppress  the  nucleation  and  growth  of  dislocation 
loops  and  voids  during  irradiation.”  LiTa03  belongs  to  a 
family  of  oxides  that  possess  a  corundum-derivative  crystal 
structure’1’  that,  because  of  their  compositional  complexity, 
are  considered  to  be  radiation  tolerant  materials.  LiTaC>3  is 
the  least  radiation  tolerant  in  this  family  of  oxides  and  is  one 
reason  it  was  chosen  for  this  work.3’ 

The  aim  of  this  paper  is  two-fold.  The  first  is  to  under¬ 
stand  the  deformation  micro-mechanisms  of  C-plane  LiTaC>3 
single  crystals  loaded  along  [0001].  The  second  is  to  explore 
the  effects  of  He  ion  irradiation  on  the  deformation 
mechanisms. 


II.  EXPERIMENTAL  DETAILS 

Two  high  quality,  (0001)  or  C-plane  orientation,  LiTa03 
single  crystals  were  purchased  (Yamaju  Ceramics  Co., 
Aichi,  Japan)  with  both  sides  polished  to  a  mirror  finish. 

The  NI  experiments  were  performed  at  room  temperature 
with  a  nanoindenter  (XP  system,  MTS  Corp.,  Oak  Ridge,  TN) 
equipped  with  a  continuous  stiffness  measurement  (CSM) 
attachment.  Three  diamond  hemispherical  indenters  with  radii, 
R,  of  21  /tm,  5  /tm,  and  1.4  //rn  were  used.  Typically,  a  tip 
was  repeatedly  indented  in  the  same  location  to  a  given  load 
along  [0001],  The  loading  rate/load  ratio  was  constant  at  0.1. 
To  correct  for  instrumental  drift,  the  unloading  segments  of 
the  sixth  and  subsequent  cycles  were  shifted  so  as  to  align 
them  with  the  corresponding  unloading  segment  of  the  previ¬ 
ous  cycles  before  the  results  were  converted  to  NI  stress-strain 
curves.  This  was  carried  out  if  and  only  successive  load-dis¬ 
placement  cycles  had  identical  areas  (see  below). 

The  load-displacement  results  were  zero-point  cor¬ 
rected.  To  determine  the  effective  zero  point,  we  used  the 
method  of  Moseson  et  al.24  The  latter  exploits  the  fact  that, 
for  a  spherical  tip,  the  following  relationship: 

S  =  2  E*a,  (7) 


where  S  and  a  are  harmonic  contact  stiffness  of  the  surface  and 
contact  radius,  respectively,  holds.  More  details  can  be  found 
in  Ref.  24.  The  effective  modulus  of  the  surface,  E*  is  given  by 


1  _  (1  -v2)  (1  -  0.072) 

E*~  E  +  1140 


(8) 


where  v  and  E  are  the  Poisson’s  ratio  and  Young’s  moduli  of 
the  sample,  respectively.  The  other  numbers  are  the  corre¬ 
sponding  values  for  the  diamond  indenter  tip.  Poisson’s  ratio 
of  LiTaC>3  is  assumed  to  be  0.25. 37 

Post-indentation  surface  features  were  examined  using 
a  scanning  electron  microscope  (SEM;  Zeiss  Supra  50VP, 
Germany).  The  Vickers  microhardness  was  measured  using 
a  load  of  10  N.  The  moduli  and  hardness  values  were  also 
measured  using  a  Berkovich  indenter  and  the  Oliver  and 
Pharr  method.20 

One  of  the  LiTa03  crystals  was  irradiated  with  2  MeV 
He4  ions  at  room  temperature  to  an  ion  fluence  of  1.67  x  1015 
He/cm2.  Ion  range,  nuclear  and  electronic  energy  partitioning 
were  estimated  using  Lindhard-Scharff-Schiott  (LSS)  proce¬ 
dure  for  calculating  ion  stopping.33^41  Based  on  LSS,  the 
range  of  2  MeV  He  1  ions  in  LiTa03  was  estimated  to  be  4.81 
/mi  (assuming  a  mass  density  for  LiTa03  of  7.41  g/cm3). 

The  lattice  damage  was  measured  by  Rutherford  back- 
scattering  spectroscopy  (RBS)  in  channeling  mode  (RBS/C) 
using  a  2  MeV  He  1  ion  beam  with  a  backscattering  angle  of 
167°.  The  ion  irradiation  and  RBS/C  analysis  were  carried 
out  at  the  Ion  Beam  Materials  Laboratory  (IBML)  at  Los 
Alamos  National  Laboratory. 

III.  RESULTS 

A.  Nanoindentation  results 

Typical  NI  load-displacement  results  obtained  when  the 
2 1  /mi  radius  tip  is  indented  into  the  unirradiated  sample  are 
shown  in  Fig.  2.  The  corresponding  results  for  the  5  jtm  tip 
radius  indenter  are  shown  in  the  bottom  right  inset  of  the 
same  figure.  The  results  for  the  irradiated  samples  were  iden¬ 
tical  to  the  unirradiated  samples  and  are  not  shown. 

In  all  load-displacement  plots  (Fig.  2),  the  first  cycle 
was  open.  After  the  indenter  was  unloaded  and  reloaded  to 
the  same  maximum  load  and  into  the  same  location,  the 
repeat  cycles  close  and  ultimately  become  fully  reversible 


FIG.  2.  (Color  online)  Typical  NI  load-displacement  curves  when  an  unirra¬ 
diated  C-plane  LiTa03  single  crystal  is  loaded  20  times  to  550  mN  with  the 
21  /.an  radius  indenter.  Note  presence  of  small  pop-ins.  Bottom  right  inset 
shows  the  same  plot  for  the  5  /. im  radius  tip  loaded  twenty  times  to  100  mN. 
In  both  cases,  for  clarity’s  sake,  only  a  few  cycles  are  plotted. 
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FIG.  3.  (Color  online)  Load-displacement  cycles  2,  3,  6,  10,  15,  and  20 
obtained  when  the  21  /an  tip  was  indented  along  the  [0001]  into  an  irradi¬ 
ated  surface.  The  curves  were  shifted  to  the  right  from  their  original  position 
for  clarity.  Cycles  2  and  3  are  open;  cycles  6  to  20  are  closed  and  equal  in 
area.  Re -loading  to  a  lower  load  after  loading  to  the  maximal  load  always 
results  in  closed,  reversible  nested  loops,  shown  for  cycle  2  only.  Inset  plots 
the  corresponding  load-displacement  loops’  areas  vs  cycle  number  obtained 
with  the  21  /mi  tip  on  an  unirradiated  surface.  After  about  6  cycles,  the  areas 
are  constant. 

and  reproducible  (Fig.  3).  However,  perfect  reproducibility 
is  only  achieved  somewhere  between  cycles  5  to  10  for  all 
locations  and  tip  sizes.  The  area  of  cycles  6  to  20  are,  within 
the  resolution  limit  of  our  NI,  identical  (inset  in  Fig.  3). 

Typical  nested  loops  obtained  when  a  given  location  is 
loaded  to  a  maximum  load,  unloaded,  and  re-loaded  to  pro¬ 
gressively  higher  loads  are  shown  on  the  loop  labeled  2  in 
Fig.  3. 

When  the  S  versus  a  results  for  the  three  radii  are  plotted 
(Fig.  4(a)),  it  is  obvious  that  the  relationship  between  the  two 
is  linear.  From  least  squares  fits  of  the  lines,  and  making  use 
of  Eqs.  (7)  and  (8),  the  average  Young’s  moduli  for  the 
1.4  fxn i,  5  /tm,  and  21  /an  tips,  were  calculated  to  be  206  ±  4 
GPa,  205  ±  5  GPa,  and  220  ±  5  GPa,  respectively.  The  for¬ 
mer  two  values  are  in  excellent  agreement  with  the  moduli 
obtained  using  a  Berkovich  tip  on  the  same  surface,  viz., 
207  ±  3  GPa. 

Why  E  for  the  21  /tm  tip  is  different  is  unclear  at  this 
time.  One  possibility  for  the  discrepancy  is  that  the  21  /an  in- 
denter  tip  may  not  be  perfectly  spherical.  In  the  case  of  the 
5  /an  and  1.4  /an  indenters,  the  total  penetration  depth  was 
kept  below  the  spherical  limit,  as  reported  by  Albayrak 
et  air1  Regardless  of  the  reasons  for  these  discrepancies,  they 
have  little  bearing  on  the  conclusions  reached  in  this  work 
that,  as  shown  below,  rely  much  more  on  the  results  obtained 
in  the  plastic  and/or  nonlinear  elastic  regimes  for  which  what 
occurs  in  the  elastic  regime  has  little  influence,  since  for  all 
intents  and  purposes  in  the  plastic  regime,  ht  ss  hc,  where  ht  is 
the  total  indentation  depth  and  hc  is  the  contact  height. 

Figure  4(b)  compares  the  stress-strain  curves  for  the 
three  tips.  In  all  cases,  the  first  cycle  delineates  two  regimes: 
a  linear  elastic  regime  followed  by  a  plastic  regime.  For  the 
1.4  /tm  and  5  /un  tips,  the  slope  of  the  elastic  regime  (shown 
by  a  dashed  inclined  line),  is  consistent  with  the  results  of  S 


a  (nm) 


Strain(a/R) 


FIG.  4.  (Color  online)  (a)  The  S  vs  a  curves  for  the  1 .4  /tm  and  5  /tm  tips, 
(b)  Typical  NI  stress-strain  curves  obtained  after  a  given  location  was 
indented  to  the  highest  loads  (550  mN  for  21  /tm,  100  mN  for  5  /tm,  and 
20  mN  for  1 .4  /tm  tips)  for  two  cycles,  unloaded  and  progressively  loaded  to 
higher  stresses  to  obtain  the  nested  loops  (three  left  curves).  Plot  shown  on 
extreme  right  shows  the  reproducible  NI  stress-strain  loops  for  1.4  /tm 
indenter;  it  was  shifted  to  the  right  from  its  original  position  for  clarity. 
Dashed  horizontal  and  inclined  lines  represent  the  Vickers  microhardness, 
Berkovich  hardness,  and  elastic  moduli  obtained  from  the  S  vs  a  curves, 
respectively.  Pop-ins  were  only  observed  when  the  21  /tm  tip  indenter  was 
used. 

versus  a  plots  (Fig.  4(a)).  As  noted  above,  the  stress-strain 
results  for  the  21  /tm  tip  indenter  in  the  elastic  regime  are 
incorrect. 

During  the  first  cycle,  in  the  plastic  regime,  the  strain¬ 
hardening  rate  is  more  or  less  constant.  More  importantly,  in 
the  plastic  regime,  the  overall  shapes  of  the  stress-strain 
curves  are  weak  functions  of  R.  Relatively,  large  observable 
pop-in  events  between  the  elastic  and  plastic  regimes  were 
only  observed  when  the  21  /mi  tip  was  used. 

At  «  6  ±  0.5  GPa,  the  Vickers  microhardness  values 
(denoted  by  a  horizontal  dashed  line  in  Fig.  4(b))  were  meas¬ 
ured  on  both  unirradiated  and  irradiated  samples  at  a  load  of 
10  N.  At  «  9  GPa,  the  Berkovich  hardness  measured  on  the 
same  surfaces  is  also  shown  as  a  horizontal  dashed  line. 

In  much  of  our  work  to  date,7’12  1  ,-|7-ISi-''-'1  q,e  defor¬ 
mation  was  initially  linear  up  to  a  pop-in  stress,  beyond 
which  the  deformation  was  plastic.  As  noted  above,  herein, 
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no  pop-ins  were  observed  for  the  1 .4  /im  tip  and  only  a  few 
for  the  5  /un  tip.  Conversely,  for  the  21  /<m  tip,  a  variation  in 
pop-in  stresses  was  observed.  The  distribution  of  pop-in 
stresses  can  be  adequately  described  by  Weibull  statistics 
(Fig.  5)  and  appear  not  to  be  a  function  of  irradiation. 

After  the  first  cycle,  two  loading  trajectories  were  fol¬ 
lowed.  The  first  was  to  load  the  same  location  20  times  to 
the  maximum  load.  Such  experiments  were  carried  out  to 
investigate  the  reproducibility  and  fully  reversible  nature  of 
the  loops  generated.  Based  on  the  results  shown  in  Fig.  4(b) 
and  the  loops  shown  on  the  extreme  right,  we  conclude  that 
these  loops  are  highly  reproducible  and  reversible. 

The  second  protocol  was  to  load  the  indenter  to  the 
highest  load  for  two  cycles,  unload,  and  then  reload  to  pro¬ 
gressively  higher  loads  to  generate  the  nested  loops  observed 
in  Fig.  4(b).8’9’25’42  Flere,  as  in  previous  work  on  LiNb03,25 
the  shape  of  the  loops  is  dependent  on  R.  The  21  /tin  and  5 
/tm  indenters  yield  shorter,  wider  loops;  the  1.4  /im  tip,  on 
the  other  hand,  results  in  elongated  thinner  loops  (Fig.  4(b)). 

All  KNE  solids  can  be  characterized  by  three  parame¬ 
ters,  cr,  eNL,  and  Wd,  all  obtainable  from  the  hysteretic  stress- 
strain  curves.  At  every  stress,  cr,  Wd,  and  £NL  were  estimated 
from  the  nested  loops,  (for  how  £NL  is  defined  in  this  work, 
see  Fig.  1(c)).  According  to  Eqs.  (4)  and  (5),  Wd  should  scale 
with  both  cr2  and  £NL,  as  observed  in  Figs.  6(a)  and  6(b), 
respectively.  The  lowest  correlation  coefficient,  6, 2  obtained 
from  least  square  analysis  of  the  results-shown  in  Fig.  6  by 
the  solid  inclined  lines-is  0.97;  most  are  >  0.98.  It  is  obvious 
that  the  model  predictions  are  well  adhered  to.  The  next  step 
is  to  try  and  quantify  some  of  the  parameters  to  ensure  that 
they  are  physically  tenable.  To  do  so,  the  following  assump¬ 
tions  were  made:  G  =  c44  =  95  GPa,43  v  =  0.25, 37  y  =  0.05, 
w  =  b  =  0.515  nm,44  kj=2,  and  M  =  2.  As  noted  above, 
experimentally  for  polycrystalline  Mg,  k,  was  found  to  vary 
between  1  and  2. 30,45,46  For  Co,  k ,  was  estimated  to  be  2. 
Reed-FIill  et  al?2  assumed  kj  =  2  for  (1121)  twins  in  Zr. 
Given  that  the  (1121)  twin  is  a  special  case  of  a  kink  bound¬ 
ary,  for  which  a  dislocation  loop  is  nucleated  every  c-lattice 


FIG.  5.  (Color  online)  Plot  of  Weibull  probabilities  (SP)  vs  pop-in  stresses 
(<r)  for  21  /im  indenter  for  both  unirradiated  and  irradiated  samples.  The 
Weibull  moduli,  m,  are  shown  on  the  figure. 


FIG.  6.  (Color  online)  Plots  of  (a)  Wd  vs  a1  and  (b)  Wd  vs  r:N,  as  a  function 
of  indenter  radius.  The  slope  dependence  on  the  tip  size  is  clear  in  (b).  Each 
line  represents  a  different  location. 

parameter,  i.e.,  D  =  c,  it  is  reasonable  to  assume  this  value 
here  as  well.47  This  comment  notwithstanding,  the  objective 
of  this  numerical  exercise  is  not  to  obtain  absolute  and  accu¬ 
rate  values  for,  say,  the  CRSS’s.  The  purpose  is  more  to 
show  that  our  model  is  consistent  with  the  results  obtained 
and  can,  in  principle,  explain  them  without  resulting  in  val¬ 
ues  that  are  not  physical.  This  is  especially  true  here,  given 
the  very  complex,  non-uniform,  state  of  stress  under  the 
nanoindenter  and  all  the  simplifying  assumptions  made. 

According  to  our  model,  the  x-axis  intercepts  of  Fig. 
6(a)  represent  the  threshold  stresses,  cr,,  needed  to  nucleate 
the  IKBs.  Using  these  cr,  values  and  the  assumptions  made 
above,  the  lengths  of  the  IKBs,  or  domain  sizes,  2a,  can  be 
estimated  from  Eq.  (1).  Once  2a  is  known,  2/lx  and  2/]y,  at 
any  cr,  can  also  be  calculated  from  Eq.  (2). 

According  to  Eq.  5,  the  slopes  of  the  lines  in  Fig.  6(b) 
should  be  equal  to  3k,Q/b.  Assuming  k,  =  2,  Q/b  can  thus  be 
calculated.  Note  that,  according  to  Eq.  (5),  the  lines  in  Fig. 
6(b)  should  go  through  the  origin,  when,  in  fact,  they  do  not. 
The  exact  reason  for  this  state  of  affairs  is  most  probably  due 
to  the  presence  of  other  non-linear  reversible  strains  that  are 
not  due  to  IKBs.  Such  strains  were  observed  when  polycrys¬ 
talline  Co  samples  were  compressed.31  Flowever,  this  dis¬ 
crepancy  is  not  believed  to  considerably  affect  the  slopes. 
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FIG.  7.  (Color  online)  SEM  images  of  NI  mark  on  unirradiated  sample 
made  with  the  (a)  1 .4  fim  tip  loaded  to  20  mN  after  first  cycle,  (b)  5  /rm  tip 
loaded  to  100  mN  after  first  cycle,  and  21  /im  tip  loaded  to  550  mN  after  (c) 
5  cycles  and  (d)  20  cycles;  (e)  same  as  (d)  but  tilted  75°.  Note  three-fold 
symmetry  of  the  linear  surface  features  best  seen  in  (b)  and  (c).  The  top  inset 
in  (a)  is  a  schematic  of  domains  forming  in  the  twins.  Top  inset  in  (b)  shows 
three-fold  symmetry  of  twins,  which  form  in  LiNbC>3,  adapted  from  Ref.  48. 
The  features  with  very  sharp  radii  of  curvature  in  (d)  and  (e)  are  kink 
boundaries. 


The  reversible  dislocation  density  can,  in  turn,  be  calcu¬ 
lated  from  Eq.  (6).  The  calculated  values  for  crt,  2a,  Nk,  Q/b, 
2/?x,  2/iy,  and  p,  the  latter  three  at  5  GPa,  are  listed  in  Table  I. 

B.  Microstructural  observations 

Despite  the  fact  that  all  of  the  tips  used  were  spherical, 
the  indentation  imprints  in  the  SEM  clearly  exhibited  three¬ 
fold  symmetry  (Fig.  7).  The  three-fold  symmetry  is  clearest 
in  Figs.  7(b)  and  7(c).  The  widths  of  the  domains  for  the 
1.4  pm,  5  pm,  and  21  pm,  shown  in  Figs.  7(a)  to  7(d), 
respectively,  were  found  to  be  a  function  of  R. 

C.  Irradiation  damage 

Almost  all  of  the  stopping  power  for  the  2  MeV  He + 
ions  in  LiTaCH  is  attributable  to  electronic  stopping.  Nuclear 
stopping  only  plays  an  appreciable  role  near  the  ion  end-of- 
range  (at  depths  from  4. 0-4. 8  pm).  The  damage  energy, 
;y(E),  which  is  the  fraction  of  the  total  energy  (E  =  2  MeV) 
that  is  consumed  in  ballistic  damage  events  (i.e.,  kinematic 
scattering  by  atomic  nuclei  in  the  target)  is  only  0.0092  MeV 
or  0.46%  of  the  primary  ion  energy. 

The  RBS  spectra  from  the  unirradiated  and  irradiated 
samples  are  shown  in  Fig.  8  and  indicated  that  the  un-irradi- 
ated  sample  was  of  a  high  quality.  The  spectrum  did  not 
change  after  irradiation  (it  was  slightly  higher,  but  within  the 
statistical  error  of  measurements),  confirming  that  the  defect 
concentration  in  the  near  surface  region  is  small  and  below 
the  sensitivity  level.  These  RBS  spectra  correspond  to  a  near 
surface  region  ~  2  pm  thick. 

IV.  DISCUSSION 

Given  that:  a)  Doukhan  et  al.6  reported  (101 2) [1011] 
twins  as  the  only  deformation  system  in  LiTaC>3  during 

Energy  (MeV) 
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FIG.  8.  (Color  online)  RBS  spectra  from  unirradiated  random,  unirradiated 
aligned,  and  irradiated  aligned.  The  Xmm  =  l-5%  (ratio  of  aligned  to  random 
yield  just  below  the  surface  peak)  indicated  very  good  quality  of  the  unirra¬ 
diated  sample.  The  aligned  spectrum  did  not  change  after  irradiation  because 
the  concentration  of  defects  in  the  near  surface  region  is  quite  small. 
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compression  along  the  c-axis;  b)  the  twinning  system  in 
LiTa03  is  identical  to  that  of  LiNb03;6  c)  the  traces  and  steps 
shown  in  Fig.  7  are  similar  to  those  observed  in  LiNb03;25 
and  d)  the  three -fold  symmetry  of  the  indentation  marks  are 
very  similar  to  those  reported  by  Park  et  al.4&~50  in  LiNb03 
(see  inset  in  Fig.  7(b)),  it  is  reasonable  to  conclude  that  the 
traces  observed  in  Fig.  7  are  due  to  the  formation  of 
(10l2)[10ll]  twins. 

In  our  previous  paper,25  we  argued  that,  since  the  load¬ 
ing  was  along  [0001],  the  distance  between  the  lines  shown 
in  Figs.  7(a),  7(b),  7(c),  and  7(d)  is  the  width  of  the  domains, 
2/1,  and  not  their  lengths,  2a.  A  schematic  of  what  we  believe 
happens  is  shown  in  the  top  inset  in  Fig.  7(a).  There  is  no 
reason  to  believe  things  are  different  here.  Also,  shown  in 
Figs.  7(d)  and  7(e)  are  some  exceedingly  sharp  bends  that 
can  only  be  due  to  kink  boundaries,  which  is  strong  indirect 
evidence  that  kinking  must  be  operative  in  this  material. 

Based  on  the  results  shown  above  and  those  listed  in 
Table  I,  and  by  analogy  with  LiNb03,  we  conclude  that 
LiTa03  is  a  KNE  solid.  As  discussed  previously,7'12-14'16'18 
the  deformation  of  KNE  solids  under  cyclic  spherical  NI  can 
be  explained  by  invoking  the  fonnation  and  annihilation  of 
spontaneously  reversible  IKB  dislocation  loops. 

According  to  the  results  shown  in  Fig.  3,  it  is  obvious 
that,  after  the  first  cycle,  hysteretic  loops  -  whose  areas,  Wd, 
get  smaller  until  they  reach  a  steady-state  value  that  is  no 
longer  a  function  of  cycling  -  evolve.  ’  ’  ’  ’  ’  These  fully 
reversible  and  reproducible  loops  are  quite  similar  to  those 
observed  in  other  KNE  solids.  The  different  parameters,  cal¬ 
culated  from  the  nested  loops,  all  yield  reasonable  results 
(Table  I). 

The  resulting  reversible  dislocation  densities  (Table  I) 
are  comparable  to  heavily  deformed  metals28  and  are,  again, 
reasonable,  considering  they  are  calculated  at  a  stress  of 
5  GPa.25  It  is  important  to  note  that,  despite  the  differences 
in  the  shape  and  size  of  the  loops  shown  in  Fig.  4(b),  at  com¬ 
parable  stresses,  p  is  a  weak  function  of  R.  The  importance 
of  this  conclusion  lies  in  the  fact  that  the  crystal  responds  to 
the  applied  stress  by  forming  dislocation  loops,  whose  total 
lengths  per  unit  volume  are  more  a  function  of  er  than  R.25 
Note,  this  conclusion  is  valid  for  the  three  tip  sizes,  even 
though  the  2/1  values  for  the  5  pm  indents  do  not  match  with 
those  of  the  1.4  and  21  pm  tips. 

The  cr,  values  shown  in  Table  I  decrease  with  increasing 
R.  As  postulated  in  our  previous  work,25  this  is  most  prob- 


TABLE  I.  Summary  of  various  measured  and  calculated  parameters  as  a 
function  of  R.  The  following  was  assumed:  y  =  0.05,  w  =  b  =  0.515  nm, 
G  =  96.8  GPa,  v  =  0.25,  M  =  2,  and  kj  -  2. 


Tip  radius  (pm) 

1.4 

5 

21 

ct(GPa)  (from  Fig.  6(a)) 

3.2  ±0.1 

2  ±  0.2 

1.8  ±0.2 

2a  (nm)  (calculated) 

114  ±5 

292  ±112 

361  ±  60 

Nkoc3 

0.45 

0.78 

1.78 

Nk(m-3) 

2.2  x  1021 

2.5  x  1020 

3.0  x  1020 

Q/b  (MPa) 

1.3  x  103 

1.1  x  103 

0.7  x  103 

2/?x  (nm)  at  5  GPa 

44  ±2 

113  ±40 

140  ±  22 

2py  (nm)  at  5  GPa 

59  ±6 

151  ±57 

1 86  ±31 

2py  (nm)  (Fig.  7.  SEM  image) 

49  ±5 

90  ±25 

180  ±40 

p  (m-2)  at  5  GPa 

4  x  1015 

3  x  1015 

5  x  1016 

ably  due  to  the  shrinking  of  domain  size,  2a.  For  smaller  R 
values,  the  domains  are  smaller  and  thus  -  from  Eq.  (1)  -  the 
threshold  stresses  should  be  larger,  as  observed  (Fig.  6  and 
Table  I). 

When  the  NI  stress-strain  responses  of  the  irradiated  and 
un-irradiated  surfaces  are  superimposed  (not  shown),  they 
were,  within  the  resolution  of  our  experiments  and  experi¬ 
mental  scatter,  identical.  Said  otherwise,  the  irradiation  did 
not  affect  the  stress-strain  curves  either  before,  or  after  the 
yield  points. 

Further  evidence  that  the  irradiation  had  little  effect  on 
the  response  is  shown  in  Fig.  5.  The  mean  and  standard  devia¬ 
tion  of  the  pop-in  stresses  for  the  irradiated  sample  was 
2.059  ±0.14;  that  for  the  unirradiated  sample  was 
2.205  ±  0.09.  Clearly,  these  values  with  their  uncertainties 
overlap.  It  follows  that  the  difference  we  see  in  Fig.  5  cannot 
be  significant.  This  is  especially  true  since  taking  these  values 
at  face  value  implies  that  irradiation  increased  the  pop-in 
stresses,  a  conclusion  which  is  difficult  to  reconcile  with  our 
previous  work  on  mica  that  showed  that  the  pop-in  stresses 
decreased  with  defect  concentration.17  The  same  is  true  here; 
it  is  difficult  to  argue  for  a  mechanism  where  irradiation  would 
somehow  suppress  the  pop- in  stresses. 

The  reason  for  this  state  of  affairs  is  believed  to  be  the 
mild  nature  of  the  irradiation.  Using  a  modified  Kinchin- 
Pease  equation  to  estimate  the  number  of  displaced  atoms 
per  incident  ion.51  we  find  that  the  number  of  displaced 
atoms  per  ion  is  ~  92  (assuming  the  displacement  threshold 
energy,  Ed,  for  all  target  atoms  -  Li,  Ta,  and  O  -  is  given  by 
Ed  =  40  e V,  which  is  an  arbitrary  assumption  often  used  as  a 
reasonable  guess  for  ceramics  when  there  are  no  measured 
values  available).  The  peak  displacement  damage  dose  is 
approximately  0.02  dpa  at  a  target  depth  of  ~  4.6  pm.  How¬ 
ever,  the  average  displacement  damage  dose  is  only  about 
0.001  dpa  over  the  first  2  pm  of  target  depth. 

In  summary,  over  the  first  two  microns  of  target  depth, 
the  2  MeV  He  1  ions  implanted  to  a  fluence,  0=  1.67  x  1015 
He/cm2,  produce  only  about  one  Frenkel  pair  per  1000  target 
atoms.  About  99.9%  of  the  ion  stopping  over  this  target  vol¬ 
ume  is  due  to  electronic  stopping.  Unless  LiTa03  is  highly 
susceptible  to  radiolysis  (i.e.,  permanent  atomic  displace¬ 
ments  due  to  electronic  excitations  and  bond  breaking),  then 
these  irradiation  conditions  will  produce  little  in  the  way  of 
permanent  point  defects  or  atomic  disorder.  This  would 
explain  why  no  detectable  changes  in  mechanical  properties 
were  observed  in  the  NI  experiments.  In  future  experiments, 
ion  fluences  that  will  produce  more  ballistic  damage  over  the 
first  two  microns  will  be  used. 

For  reasons  that  are  not  entirely  clear,  twins  are  easier  to 
nucleate  when  the  sharper  tips  are  used.  This  is  best  seen  by 
the  lack  of  pop-in  stresses  when  the  1.4  and  5  pm  tips  are 
used.  With  the  21  pm  tip,  the  strain  energy  has  to  reach  a  cer¬ 
tain  value  before  the  twins  are  nucleated.  This  is  an  impor¬ 
tant  observation  that  needs  to  be  further  looked  into,  since  it 
signifies  that  the  nucleation  of  twins  depends  on  more  than 
the  stress  under  the  indenter. 

In  our  work  on  the  MAX  phases,  we  have  shown  that 
Q/b  or  the  CRSS  of  the  IKB  dislocations  is  inversely  pro¬ 
portional  to  the  square  root  of  the  grain  or  domain  size.10 
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To  check  whether  this  is  true,  here  we  plotted  Q/b  versus 
\/j2<x  (Fig.  9).  A  decent  correlation  was  found.  If  one 
plots  Q/b  versus  1/2)3,  measured  directly  from  the  SEM 
micrographs  (see  below),  an  even  better  correlation  is 
found  (Fig.  9).  Since  a  is  proportional  to  /?,  this  correlation 
again  is  consistent  with  Q/b  values  that  follow  a  Hall- 
Petch  like  relationship.  Lastly,  one  can  also  plot  Q/b  ver¬ 
sus  1/yR;  again  the  correlation  is  acceptable.  Note  that 
the  correlation  found  in  Fig.  9  is  not  a  result  of  a  circular 
argument,  since  a  is  calculated  from  ax  and  the  CRSSs  are 
calculated  from  the  totally  independent  measurements  of 
Wd  and  sNL.  The  correlations  between  Q/b  versus  l/x  R  or 
1/^2/?  are,  thus,  also  totally  independent. 

When  Q/b  is  plotted  versus  1A/R  for  the  LiNb03  single 
crystals,  -  shown  as  a  black  dashed  line  in  Fig.  9  -  using 
the  methodology  used  herein,  the  correlation  is  once  again 
excellent.  (The  calculations  of  Q/b  made  in  Ref.  25 
made  slightly  different  assumptions  to  reach  different  values 
of  Q/b).  Not  only  is  the  correlation  excellent,  but  the  result¬ 
ing  line  is  almost  parallel  to  that  for  LiTa03  (the  line  at 
extreme  left  in  Fig.  9).  This  is  an  important  result  for  the 
following  reason:  if  one  makes  the  reasonable  assumption 
that  the  CRSSs  of  the  IKB  dislocations  are  proportional  to 
c44,  then  the  ratio  of  CRSSs  for  LiTa03  and  LiNb03  should 
be  rs  1.62.  Gratifyingly,  the  ratio  of  CRSSs  -  for  a  given 
domain  size  denoted  by  the  vertical  dashed  line  -  obtained 
from  the  results  shown  in  Fig.  9  is  ss  1.6.  Whether  this  is 
coincidental  or  not  needs  more  work  on  other  crystals  with 
different  c44  values,  but  is  certainly  consistent  with  the  ideas 
proposed  herein. 

Comparing  both  calculated  and  measured  values  for  2)3y 
in  Table  I  shows  that  the  calculated  values  from  the  model 
and  measured  values  from  SEM  images  for  1.4  /mi  and 
21  /mi  are  in  good  agreement  with  each  other.  For  reasons 
that  are  unclear,  for  the  5  /mi  indents,  the  calculated  values 
differ  from  the  measured  ones.  Also,  the  calculated  length  of 
the  IKBs,  2a,  for  1.4  /mi,  5  /mi,  and  21  /mi  indents  are  quite 
reasonable  and  scale  with  the  indenter  radii. 


FIG.  9.  (Color  online)  Hall-Petch-like  correlation  between  the  CRSS  of  the 
IKB  dislocations  (fl/b)  and  1  /v  domain  size,  where  the  domain  size  is 
assumed  to  be  2a,  2/?,  or  R.  The  black  dashed  inclined  line  represents  Q/b  vs 
1/^/R  for  LiNb03  (see  Ref.  25). 


Based  on  the  totality  of  our  results,  the  following  sce¬ 
nario  for  what  happens  under  the  spherical  NI  tips  can  be  re¬ 
created.  At  pop-in  for  the  21  /tm  indenter,  or  at  the  yield 
points  for  the  other  two  tips,  twins  form.  These  twins,  in 
turn,  rotate  basal  planes  into  orientations  that  are  more  ame¬ 
nable  for  basal  slip.  Concomitantly,  or  as  a  result,  the  single 
crystal  is  fragmented  into  much  smaller  domains,  with  a  size 
that  scales  with  R.  At  a  threshold  stress  -  that  is  inversely 
proportional  to  the  square  root  of  the  domain  size  -  IKBs  nu¬ 
cleate  within  these  domains.  The  IKBs  are  fully  reversible, 
and  the  to-and-fro  motion  of  the  IKB  dislocations  dissipates 
substantial  amounts  of  energy. 

Note  that  IKBs  cannot  form  in  single  crystals  because, 
once  they  become  critical,  they  would  simply  run  to  the  ends 
of  the  crystal  and  devolve  into  mobile  dislocation  walls.11 
Equation  (1)  -  derived  using  a  Griffith-like  approach  -  is 
based  on  that  scenario.  IKBs  can  thus  only  form  when  there 
are  grain,  twin,  or  domain  boundaries  that  confine  them.  Dur¬ 
ing  NI,  these  domains  form  during  the  first  cycle,  especially 
after  pop-in  events  or  yield  points.  It  was  expected  that  irra¬ 
diation  defects  on  the  surface  reduce  the  pop-in  stress. 

We  note  that  the  Vickers  hardness  values  and  the  min¬ 
ima  in  the  stress-strain  curves  after  pop-ins  or  simply  the 
yield  points,  if  pop-ins  are  not  present  (Fig.  4(b)),  are  in 
excellent  agreement.  This  is  an  important  result  since  it 
indirectly  confirms  our  methodology  for  converting  NI  load/ 
displacement  to  NI  stress-stain  curves.  The  same  observation 
was  made  in  most  of  our  other  previous  work.12'14'21’  The 
fact  that  the  Berkovich  modulus  is  also  in  good  agreement 
with  our  value  is  another  independent  confirmation  for  our 
NI  methodology. 

Lastly,  it  is  worth  emphasizing  that,  given  the  many 
simplifying  assumptions  made  in  our  KNE  model,  such  as 
assuming  a  uniform  uniaxial  stress  state,  which  is  far  from 
what  is  happening  during  NI,  to  the  definition  of  strain  to  be 
a/R  and  other  simplifying  assumptions,  the  agreement 
between  theory  and  experiment  has  to  be  considered  excel¬ 
lent.  One  reason  for  this  state  of  affairs  is  that  the  calculation 
of  contact  stresses  is  reasonably  straightforward.  Moreover, 
since  the  non-linear  strains  are  small,  the  results  obtained  are 
a  weak  function  of  the  exact  definition  of  strain.  These  com¬ 
ments  notwithstanding,  it  is  hereby  acknowledged,  as  noted 
above,  that  the  absolute  values  of  CRSS  calculated  herein 
have  to  be  taken  with  a  large  grain  of  salt  because  of  all  the 
aforementioned  simplifications,  etc. 

V.  CONCLUSIONS 

The  deformation  response  of  C-plane  LiTa03  single 
crystals  was  studied  using  spherical  indenters  with  three  dif¬ 
ferent  radii.  When  the  load/displacement  results  are  con¬ 
verted  to  NI  stress-strain  curves,  we  conclude  that: 

(a)  During  the  first  NI  cycle,  two  regimes  are  observed:  a 
linear  elastic  regime  followed  by  a  plastic  regime,  in 
which  strain-hardening  is  observed.  For  the  smaller  tips, 
plastic  deformation  is  accompanied  by  (1012)  twins.  In 
the  case  of  the  21  /mi  indenter,  the  elastic  regime  was 
separated  from  the  plastic  regime  by  pop-ins.  At  the  pop- 
ins,  twins  are  nucleated. 
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(b)  The  elastic  moduli  determined  from  S  versus  a  curves  of 
the  1.4  /mi  and  5  /tm  tips  was  205  ±  5GPa.  The  corre¬ 
sponding  Berkovich  modulus  is  207  ±  3  GPa. 

(c)  No  noticeable  effect  of  the  irradiation  was  observed  on 
the  NI  results  of  the  sample.  The  main  reason  for  this 
state  of  affairs  is  believed  to  be  the  mildness  of  the  irra¬ 
diation  conditions. 

(d)  When  spherical  indenters  were  repeatedly  loaded  to  a 
given  maximum  load,  fully  reversible,  reproducible  hys¬ 
teresis  loops  are  obtained.  The  G/b  values  calculated 
from  the  model  were  inversely  proportional  to  the  square 
root  of  the  domain  size. 

(e)  The  ratio  of  the  CRSSs  for  LiNb03  and  LiTaOa  deter¬ 
mined  from  our  IKB -based  microscale  model  were  found 
to  be  roughly  equal  to  the  ratios  of  the  shear  moduli  of 
these  two  compounds. 
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